Theory of strong localization effects of light in disordered loss or gain media 
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We present a systematical theory for the interplay of strong localization effects and absorption or 
gain of classical waves in 3-dimensional, disordered dielectrics. The theory is based on the selfcon- 
sistent Cooperon resummation, implementing the effects of energy conservation and its absorptive 
or emissive corrections by an exact, generalized Ward identity. Substantial renormalizations are 
found, depending on whether the absorption/gain occurs in the scatterers or in the background 
medium. We find a finite, gain-induced correlation volume which may be significantly smaller than 
the scale set by the scattering mean free path, even if there are no truly localized modes. Possible 
consequences for coherent feedback in random lasers as well as the possibility of oscillatory in time 
behavior induced by sufficiently strong gain are discussed. 
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I. INTRODUCTION 

The strong or Anderson localization (AL)i of light in 
d = 3 dimensional, random media has remained a fasci- 
nating and controversial issue, despite its early theoreti- 
cal anticipation-ifii and despite the experimental observa- 
tion of its precursor effect, the coherent backscattering 
light cone of weak localizationi4 Early experimental re- 
ports of an anomalously low diffusion constant in strongly 
scattering powders^ were explained by a reduction of the 
transport velocity by Mie resonances A However, recent 
works on the experimental proof of light localization^ has 
lead to a controversial, still ongoing debate The reason 
resides in the difficulty to distinguish AL from absorption 
in the medium. The problem becomes even more pressing 
for localization in random media with stimulated emis- 
sion {random lasers) Here the scattering mean free 
path is large compared to the wave length, i X, so 
that localization effects are expected to be small. How- 
ever, the experiments indicate unambiguously, by direct 
observationii and especially by measuring the photon 
statisticsfi^ that the laser emission is due to coherent 
feedback and occurs from spatially confined spots in the 
sample. As an alternative to lasing from Anderson lo- 
calized modes, the existence of preformed cavities in the 
random medium has been proposed in Ref. Il3l 

AL has been understoodi^ii^ as an effect of repeated 
self-interference (so-called Cooperon contributions) of 
diffusive modes. Since diffusion, as a hydrodynamic phe- 
nomenon, relies on particle number or energy conserva- 
tion, this raises the fundamental question for the fate of 
AL in active media altogether. Intuitively one expects 
coherent amplification to enhance the transmission, but 
at the same time also the Cooperon interference. Since 
the latter tends to localize a wave and hence to reduce 
the transmission, the overall effect is left an open issue. 
Paasschens et al.^^ , treating the transmission in a lin- 
earized way, prove analytically in c? = 1, that amplifi- 
cation diminishes the transmission just like absorption 
does. In contrast, Jiang et alJ^, going beyond the lin- 



ear approximation, find numerically an enhanced trans- 
mission in the long-time limit. Although important in- 
sight has been obtained by further numerical studies in 
d — l^^ and d — 2,^^ this controversial situation calls for 
a semi-analytical theory to systematically analyze trans- 
port quantities like the diffusion coefficient D and the 
intensity correlation length ^a- This is especially so in 
d = 3 where the realization of AL modes is unclear even 
in the passive system. 

In this article we present such a theory, based on the 
selfconsistent Cooperon resummation pioneered by VoU- 
hardt and Wolfle.^^ It allows to distinguish AL proper- 
ties from absorption/gain-induced decay or growth. The 
evaluations will be done for d = 3, although the the- 
ory is valid for d = 1, 2, 3. As the central results we (1) 
recover, for small gain, the duality of AL with respect 
to absorption or gain, together with a spatially homoge- 
neous, exponential decay/growth in the long-time limit, 
thus reconcihng the results of Refs. ll6lf7i (2) In addition, 
we find strong loss or gain dependent renormalizations of 
D and of ^a, asymmetrical in loss or gain, originating not 
from wave interference but from the violation of the con- 
servation laws. These renormalizations are, in particular, 
sensitive to whether the loss (gain) occurs in the homoge- 
neous background medium or in the random scatterers, 
and thus do not occur in d = 1 layered structures. (3) 
We predict that due to loss or gain the intensity corre- 
lation length can be substantially smaller than the scat- 
tering mean free path £, constituting a finite coherence 
volume. Possible implications on coherent feedback in 
random lasers are discussed. 



II. MODEL AND TRANSPORT THEORY 

We consider a system of randomly positioned spherical 
scatterers in a background medium with dielectric con- 
stants 65 and £(,, respectively. On the semiclassical level, 
linear absorption and the stimulated part of emission are 
represented by a positive or negative imaginary part of 



2 



the index of refraction, = ^/el ' ~ ■ Since 
we will focus on the fundamental problems posed above 
rather than a quantitative description, we will neglect the 
vector nature (polarization) of light and consider scalar, 
classical waves only. These obey the wave equation 
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e(r)*c^(r) + V^'J'„(f) = -iuj—j^{r) 



(1) 



where c denotes the vacuum speed of light and ju{r) 
an external source. The dielectric constant e(r) = + 
Ae V{r), Ae = eg — Cb, describes the arrangement of scat- 
terers through the function ^(r ) = J2r '^r )j '^ith 

(f ) a localized shape function at random locations R. 
The Fourier transform of the retarded, disorder averaged 
Green's function of Eq. (Q) reads. 



The right-hand side of Eq. Q represents reactive effects 
(real parts), originating from the explicit w^-dependence 
of the photonic random "potential" . In conserving me- 
dia (Imcb = hnCs = 0) these terms renormalize the 
energy transport velocity relative to the average 
phase velocity Cp without destroying the diffusive long- 
time behaviorisms In presence of loss or gain, however, 
these effects are enhanced via the prefactor fui{^) = 
{u}QReAe + iuj'^lmAe)/{uj^ReAe + iu}QlmAe), which now 
does not vanish in the limit ^ 0. As a consequence, 
more severe renormalizations of the diffusion coefficient 
D and a finite intensity correlation length S,a are induced, 
as seen in section Unl 



III. SELFCONSISTENT SOLUTION 



" efc(c./c)2 - |fc|2 



(2) 



where the retarded selfenergy E| arises from scattering 
off the random "potential" -(w/c)2(e, - ib)V{f). The 
mode density is given by N{uj) = — (aj/7r)ImG'o , Gq = 
Jd'k/{27rrGl 

To develop a transport theory we now turn to the 4- 
point intensity correlation function, defined in terms of 
the non-averaged Green's functions G, G* and the dis- 
order average (...) as $^,-,(Q,f^) = (G^;,-; ^ ) . It is 
determined by the kinetic equation, see, e.g., Ref. l20l 
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Here we have introduced the usualSSi center-of-mass and 
relative frequencies and momenta: The variables Q, Q 
are associated with the time and position dependence 
of the averaged energy density, with Q = Q/\Q\, while 
Lu± = w ± r2/2 and q± = qzLQ/2 etc. are the frequencies 
and momenta of in- and out-going waves, respectively. 
In order to analyze the correlation function's long-time 
(ri 0) and long-distance {\Q\ 0) behavior, terms of 
0{fl^, Q^, ^Q) have been neglected here and throughout 
this paper. Eq. Q contains both, the total quadratic mo- 
mentum relaxation rate l/r^ = Im(ebw2/c2 — E") (due 
to absorption/gain in the background medium as well as 
impurity scattering) and the irreducible two-particle ver- 
tex function 7~, ((5,r2). The energy conservation is im- 
plemented in a field theoretical sense by a Ward identity 
(WI) which has been derived for the photonic case in Ref. 
I20I and which for scalar waves takes the exact form 
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As long as wave interference effects are not taken into 
account, the selfenergy and the irreducible vertex may 
be evaluated within a local approximation, like the in- 
dependent scatterer approximation employed in Ref. |^ 
or modifications of the coherent potential approxima- 
tion (CPA). 22 In this case, these quantities become q- 
independent, ^ E^, 7^, .(Q, rj) ^ 7^, and Eq. @ 
may be solved in a straight-forward way^*? with the help 
of Eq. However, to account for Anderson localiza- 

tion effects, the ^-dependence of 7^, -(Q, ^i) is known to 
be essentialii^ii^ To solve the kinetic equation Q for this 
case, <I>^^ is expandedSS into its moments with respect to 

the longitudinal current vertex {cpp- Q), i.e. into the en- 
ergy density correlator P^{Q, ri) and the energy current 
density correlator J^((3, ri),^-^" 
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In these definitions the phase velocity Cp of the disordered 
medium is defined via the dispersion obtained from the 
averaged single-particle Green's function, Eq. as 



Cp = Re 



(7) 



and the transport velocity We will be given below. In- 
serting the expansion by moments into the kinetic equa- 
tion ||2J) and employing the WI, one obtains the exact 
continuity equation and the current relaxation equation, 
respectively, 
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J^+AQP^ = 0, (9) 
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The expression for the transport velocity is now deter- 
mined by the condition that the density and the current 
correlators must appear on the left-hand side of the con- 
tinuity equation ((HJ without additonal prefactors. One 
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The expressions for the coefficients A, A, A, g, 
follow in analogy to Ref . |23, 
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They enter into the transport quantities given explicitly 
below. Eq. ^ represents the local energy conservation 
with corrections due to absorption or gain as well as re- 
active corrections. Eq. describes the time evolution 
of the current density ( ) induced by a density distribu- 
tion (P^) due to momentum relaxation and the memory 
kernel M{Vi). M{Vi) is a functional of the irreducible 2- 
wave vertex 7|^/^(Q,S^)- Its characteristic f2-dependence 
describes memory effects due to the diffusion and inter- 
ference as well as amplification or absorption of waves in 
the medium, see below. Eqs. © are easily solved to 
obtain the diffusion pole form for the density correlator, 
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with the il-dependent diffusion coefficient -D(ri), 

D{n) [I - i rJcjT^Rceb] = Dl°^ - c^T^D{n)M{Lj) (12) 

and the absorption or gain induced correlation length 
of the diffusive modes. 
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The diffusion constant without memory effects, Dq°* - 
Dh + Z?s, consists of the bare diffusion constant,'^ 
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and renormalizations from absorption or gain in the back- 
ground medium {Di,) and in the scatterers {Dg), 



= lr,A,T^Do, 



(15) 
(16) 
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FIG. 1: (Color online) Normalized diffusion constant D{0) for 
£(, = 1 and scatterer volume fraction u — 50%. The frequency 
unit in all figures is uiq ~ 27rc/ro with ro the scatterer radius. 
The localization transition is reached for es > 24. Inset: phase 
diagram in the plane within the first Mie resonance for 

£(, = 1 and various e^. Shaded areas: localized phases. 



where Dq is the same as in Eq. (|14f) . with (ImG^)^ re- 
placed by (ReG| 2). In Eqs. (HSJl-llinil the following short- 
hand notations have been introduced, 



A, 

Be 



Im(AeS") 
Im(AeG^) ' 
2[u,ReGo + ReI]o] 
(ReAe)2 + (ImAe)^ 
2w2(ReAe)2 ' 



ImAe/ReAe, 



We have carefully analyzed 7^', - for the selfconsistent 
calculation of M(ri) fi^424 exploiting time reversal sym- 
metry of propagation in the active medium. In the long- 
time limit (f2 0) the dominant contributions to 7^/g- 
are the same maximally crossed diagrams (Cooperons) as 
for conserving media, however now acquiring the absorp- 
tion (gain)-induced decay (growth) rate ^a^D. M{il) 
reads 
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Eqs. lfT^ - lfT7|) constitute the selfconsistency equations for 
the diffusion coefficient D{n) and the correlation length 
in presence of absorption or gain. We have evaluated 
them below em plo ying the independent scatterer approx- 
imation of Ref. I21I for the single-particle quantities. More 
sophisticated approximation schemes, suitable for high 
scatterer concentrations, like the coherent potential ap- 
proximation (CPA), 22 will be considered elsewhere. 
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FIG. 2: (Color online) (a) Diffusion constant D{0) for u = 
50%. The scatterers exhibit different amounts of absorption 
as indicated, (b) Absorption length in units of the scatterer 
radius for the systems in (a). 
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FIG. 3: (Color online) Diffusion constant D{0) for absorbing 
scatterers {u — 50%) in a host medium with optical gain, as 
indicated. Parameters are chosen such that > 0. 



IV. DISCUSSION 

Let us first discuss thie kernel Af(il). It describes 
the enhanced backscattering. For a conserving medium 
(Imefc — Imes = 0), where — 0, it drives the 

AL transition due to its negative, infrared divergent 
contribution. '^'^^ For illustration we show in Fig. the 
AL phase diagram and D{Q = 0) for = 1 and various 
real values of Cs , displaying strong suppression near the 
Mie resonances, albeit for rather high dielectric contrast 
Ae. Introducing now absorption or gain, the quadratic 
inverse correlation length becomes non-zero and can 
assume both positive and negative values. True AL due 
to repeated enhanced backscattering is no longer possible 
in this case, because, through selfconsistency, a vanishing 
£1(0) would drive ^"^ ^ oo and D{0)M{0) 0. Inspec- 
tion of M(ri) shows analytically that to leading order in 
the real part of the integral in Eq. (|17|l is symmet- 
rical with respect to "^oT^i AL is suppressed 
by absorption or gain in a symmetrical way. This is in 
agreement with the surprising absorption/gain duality of 
which is valid for small gain and short times. 

However, we do find important deviations from this 
result for systems where there is no symmetry between 
background and scattering medium (like in our d = 3 
case): (1) As seen from Eq. (fT^ . itself is not sym- 
metric in the sign of Imeb or Ime^ and can be positive 
even for purely emissive media, Ime^ < 0, Ime^ < 0. 
(2) The full diffusion coefficient has additional contribu- 
tions Dh, Ds from loss/gain in the background and in 
the scatterers, Eqs. (|15|) . which to our knowledge have 
not been reported before. While Db is always positive 
for absorption and negative for gain as expected, Dg has 
a complicated dependence on the signs of Ime;,, Imeg and 
depends sensitively on whether absorption and/or gain 
occurs in the background medium or in the scatterers. 
This is because, in contrast to I?^, the impurity scatter- 
ing contribution Dg results from an intricate interplay be- 



tween elastic momentum relaxation and absorption/gain 
processes. We emphasize that the existence and the form 
of Df,, Ds, and are a direct consequence of the non- 
conserving terms in the WI Eq. Q), and, thus, are exact. 

The complete scenario of localization effects is now 
as follows. Even though in the presence of absorption 
or gain there are no true Anderson localized modes be- 
cause of the finite the contributions Df, and Dg 
can strongly suppress the total diffusion coefficient D{Q). 
This is shown in Fig. [3 for a system of absorbing scat- 
terers embedded in air. Even moderate absorption dras- 
tically decreases the diffusion constant D{0) close to the 
low-order Mie resonances. At the same frequencies the 
correlation length S,a is suppressed even more dramati- 
cally. For the case of a purely absorbing system, S,a may 
be identified with the effective absorption length for dif- 
fusive modes. For the case of emission, e.g. in the back- 
ground, two scenarios are possible, (i) > 0; this may 
occur due to the subtle interplay with momentum relax- 
ation processes and absorption in the scatterers. Both 
D{0) and the intensity correlation length are real and 
finite, but suppressed near the Mie resonances, as shown 
in Fig. 121 It is also seen that, e.g., absorption in the 
scatterers can be partially compensated by emission in 
the background, (ii) Re < 0; this is reahzed for 
sufficiently strong gain. It implies a pole on the real axis 
in the integration range of Eq. (|17|l and, hence, complex 
£'(0) and ^a- Fourier transforming Eq. |[TT|) . this means 
an exponential intensity growth for long times with rate 
I/tq = Ke[{^~'^+Q^)D{Q)], in qualitative agreement with 
Ref. JJ7, thus reconciling these long-time results with the 
weak gain or short-time results of Ref. This growth 
is modulated by temporal and spatial oscillations with 
characteristic frequency iln — —Q^lmD{0) and wave 
number kjj = Iml/^a, respectively. We interpret this 
oscillatory behavior as a memory effect, originating from 
the competition between the enhanced backscattering of 
waves and their amplified propagation in the surrounding 
medium. Selfinduced oscillations have been found before 
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in other driven systems with competing dynamics. 

V. CONCLUSION 

In conclusion, we have presented a semianalytical the- 
ory for the interplay of strong localization effects and ab- 
sorption or gain. True AL is not possible in the presence 
of either loss or gain. However, strong renormalizations 
of the diffusion constant D arise from the violation of the 
conservation laws. These renormalizations depend sensi- 
tively on whether absorption/gain occurs in the scatter- 
ers or the background. Intimately connected with the 
suppression of D is the appearance and reduction of a fi- 
nite intensity correlation length i even though there are 
no truly localized modes. S,a includes effects of both, im- 
purity scattering and loss/gain in the medium, and, thus, 
can be shorter than the scattering mean free path £. For 
example, in a pure medium with loss/gain would char- 
acterize the absorption/gain length and would be finite, 
while £ would obviously be infinite. It should be em- 
phasized that the present theory incorporates both, the 
physics of multiple impurity scattering valid at length 



scales larger than £ (diffuson and Cooperon contribu- 
tions) and the physics of coherent amplification present 
at all length scales. Therefore, it is expected to give an 
accurate estimate for S^a- We conjecture that in random 
lasers this finite length scale might define the coher- 
ence volume necessary for resonant feedback, that is ob- 
served experimentally^^ and that appears to be smaller 
than £ in those experiments. In order to substantiate 
this conjecture, the present localization theory should be 
coupled selfconsistently to the laser rate equations, which 
will yield results for the position dependent dielelectric 
function above the lasing threshold (Ime^.s ^ 0) and, 
hence, for the lasing mode volume. This will be sub- 
ject of further research. Oscillatory behavior in space 
and time is also predicted for sufficiently strong gain, al- 
though it is presumably difficult to observe, as it occurs 
only during the exponential intensity growth between the 
laser threshold and saturation. 
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